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Coulomb drag between two unhybridized graphene sheets separated by a dielectric spacer has 
recently attracted considerable theoretical interest. We first review, for the sake of completeness, the 
main analytical results which have been obtained by other authors. We then illustrate pedagogically 
the minimal theory of Coulomb drag between two spatially-separated two-dimensional systems of 
massless Dirac fermions which are both away from the charge-neutrality point. This relies on second- 
order perturbation theory in the screened interlayer interaction and on Boltzmann transport theory. 
In this theoretical framework and in the low-temperature limit, we demonstrate that, to leading {i.e. 
quadratic) order in temperature, the drag transresistivity is completely insensitive to the precise 
intralayer momentum-relaxation mechanism {i.e. to the functional dependence of the scattering 
time on energy). We also provide analytical results for the low-temperature drag transresistivity 
for both cases of "thick" and "thin" spacers and for arbitrary values of the dielectric constants 
of the media surrounding the two Dirac-fermion layers. Finally, we present numerical results for 
the low-temperature drag transresistivity in the case in which one of the media surrounding the 
Dirac-fermion layers has a frequency-dependent dielectric constant. We conclude by suggesting an 
experiment that can potentially allow for the observation of departures from the canonical Fermi- 
liquid quadratic-in-temperature behavior of the transresistivity. 
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I. INTRODUCTION 

Electron-electron interactions are a source of coupling 
between closely spaced nano-electronic circuits. This 
coupling has commanded a great deal of attention during 
the past thirty years or so, since it constitutes a poten- 
tial alternative to the inductive and capacitive couplings 
of conventional electronics. Early on it was realized^"^' 
that "Coulomb mutual scattering" between spatially sep- 
arated electronic systems provides a mechanism to relax 
momentum that tends to equalize drift velocities. This 
intrinsic friction due to electron-electron interactions is 
modernly referred to as "Coulomb drag'^^ffl. Early ex- 
perimental work was carried out by Gramila et al.^'^ and 
by Sivan, Solomon, and ShtrikmarpJ^ in semiconductor 
double quantum wells. 

In these experiments a constant current is imposed on 
the two-dimensional (2D) electron gas in one of the wells 
(the "active" or "drive" layer). If no current is allowed 
to flow in the other well (the "passive" layer), an elec- 
tric field develops whose associated force cancels the fric- 
tional drag force exerted by the electrons in the active 
layer on the electrons in the passive one. The transre- 
sistance pr>, defined as the ratio of the induced voltage 
in the passive layer to the applied current in the drive 
layer, directly measures the rate at which momentum 
is transferred from the current-carrying 2D electron gas 
to its neighbor. Coulomb drag is ultimately caused by 
fluctuations in the density of electrons in each layer since 
two-dimensional layers with uniformly distributed charge 
will not exert any frictional forces upon each other3-l 

The study of Coulomb-coupled 2D systems has now 
been revitalized by advances which have made it pos- 
sible to prepare robust and ambipolar 2D electron sys- 
tems (ESs), based on graphene^ layers or on the surface 
states of topological insulators (TIs)P^, that are described 



by an ultrarelativistic wave equation instead of the non- 
relativistic Schrodinger equation. 

Single- and few-layer graphene systems can be pro- 
duced, for example, by mechanical exfoliation of thin 
graphiteI31 or by thermal decomposition of silicon car- 
bidei^. Isolated graphene layers host massless-Dirac 
two-dimensional electron systems (MD2DESs) with a 
four-fold (spin x valley) flavor degeneracy, whereas 
topologically-protected MD2DESs that have no addi- 
tional spin or valley flavor labels appear automaticalljffS 
at the top and bottom surfaces of a three-dimensional 
(3D) TI thin film. The protected surface states of 3D 
TIs are associated with spin-orbit interaction driven bulk 
band inversions. 3D TIs in a slab geometry offer two 
surface states that can be far enough apart to make 
single-electron tunneling negligible, but close enough for 
Coulomb interactions between surfaces to be important. 
Unhybridized MD2DES pairs can be realized in graphene 
by separating two layers by a dielectric?^ (such as AI2O3) 
or by a few layers of a one-atom-thick insulator such 
as BNi^Isni, In both cases interlayer hybridization is 
negligible and the nearby graphene layers are, from the 
point of view of single-particle physics, isolated. Iso- 
lated graphene layers can be also found on the surface 
of bulk graphitc^^ ^^ and in "folded graphene'^ (a nat- 
ural byproduct of micromechanical exfoliation), or pre- 
pared by chemical vapor deposition^^. We use the term 
double-layer graphene (DLG) to refer to a system with 
two graphene layers that are coupled only by Coulomb 
interactions, avoiding the term bilayer graphene which 
typically refers to two adjacent graphene layers in the 
crystalline Bernal-stacking configuratiorfi^l. 

DLG and TI thin films are both described at low en- 
ergies by a Hamiltonian with two MD2DES&^^ coupled 
only by Coulomb interactions, in the absence of single- 
particle tunneling. Coulomb drag between two spatially- 



separated MD2DESs has recently attracted a great deal 
of theoretical intereslP^l'^Sl xhe calculations in Refs. [24l - 
[2^ refer to the regime in which both layers are either 
electron- or hole-doped. The Coulomb drag transresis- 
tivity in this case is negative and vanishes like T^ at low 
temperatures. Despite the considerable amount of work 
published on the subject recentljEStSl!^ no clear consen- 
sus exists on the dependence of the drag transresistivity 
in the Fermi-liquid regime on carrier densities in the two 
layers, on the interlayer distance, and on the dielectric 
constants of the media surrounding the two layers. The 
main analytical results obtained earlier by other authors 
will be summarized below in Sect. Illl 

The Coulomb drag transresistivity between two 
MD2DESs in the regime in which one layer is electron 
doped and the other is hole doped has been recently cal- 
culated by Mink et al^. In this intriguing regime, the 
authors of Ref. 30 have found that pD grows logarith- 
mically upon lowering the temperature T towards the 
critical temperature Tc for exciton condensatiorPS (con- 
densation of electron-hole pairs in a dipolar condensate). 

Coulomb drag between two graphene sheets has been 
recently measured by Kim et al!^. This first experimen- 
tal study represents an important milestone since the au- 
thors of this work have shown that the "strong-coupling" 
regime, i.e. the regime in which the interlayer distance d 
is much smaller that the typical separation between two 
electrons in each layer, is easy to achieve experimentally 
with two, independently-contacted, graphene sheetsl^. 
This study has indeed fueled the recent theoretical inves- 
tigations of Coulomb drag between two MD2DESs men- 
tioned above. 

In this Article we present in a pedagogical fashion the 
minimal theory of Coulomb drag between two spatially- 
separated MD2DESs in the regime in which both layers 
are either electron- or hole-doped. We will be only con- 
cerned with the so-called "Fermi-liquid regime" in which 
both layers are away from the charge neutrality point. 
Our theory relies on second-order perturbation theory 
in the screened interlayer interaction and on Boltzmann 
transport theory. In this theoretical framework and in 
the low-temperature limit, we demonstrate that, to lead- 
ing {i.e. quadratic) order in temperature, the drag tran- 
sresistivity is completely insensitive to the precise in- 
tralayer momentum-relaxation mechanism [i.e. to the 
functional dependence of the scattering time on energy). 
This is in contradiction with the findings reported in 
Refs. |27] and [28l We also provide new analytical results 
for the low-temperature drag transresistivity in both 
cases of "thick" and "thin" spacers, correcting in the lat- 
ter case a mistake contained in Ref. [3H1 At odds with all 
the previous literature, our results hold true for arbitrary 
values of the dielectric constants of the media surround- 
ing the two Dirac-fermion layers. Finally, we present 
numerical results for the low-temperature drag transre- 
sistivity in the case in which one of the media surround- 
ing the two MD2DEs has a strongly-frequency-dependent 
dielectric constant. We conclude by suggesting an ex- 
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FIG. 1: (Color online) A side view of the double-layer system, 
which explicitly indicates the dielectric model used in these 
calculations. The two layers hosting massless Dirac fermions 
are located &t z = Q and z = d. In a Coulomb-drag transport 
setup a constant current flow is imposed in one layer (the 
bottom one, say). If no current is allowed to flow in the other 
layer, an electric fleld develops whose associated force cancels 
the frictional drag force exerted by the electrons in the bottom 
layer on the electrons in the top one. 



periment with a DLG deposited on SrTiOgp^l that can 
pave the way for the observation of departures from the 
canonical Fermi-liquid quadratic-in-temperature behav- 
ior of the transresistivity. 

Our manuscript is organized as follows. In Sect, [ill we 
report a summary of the main analytical results which 
have been obtained by other authors. In Sect. |III| we 
present the model Hamiltonian and the most important 



basic definitions. In Sect. IV we present the Kubo formal- 
ism approach to the calculation of the drag conductivity, 
while in Sect. |V] we present a series of simplifications that 
lead to the Boltmann-transport expression for the drag 
conductivity and resistivity. These two Sections do not 
contain original results but make the paper completely 
self-contained. Experts can skip Sects. |IV||V| and go di- 
rectly to Sects. VlllVII which contain the most impor- 
tant results of this work and all our original results. In 
Sect. |VlII| we summarize our main findings and draw our 
main conclusions. 



II. SUMMARY OF THE MAIN ANALYTICAL 

RESULTS OBTAINED EARLIER BY OTHER 

AUTHORS 

Despite the large body of theoretical work ded- 
icated to Coulomb drag between two unhybridized 
MD2DESd24H29|^ no consensus appear to exist among dif- 
ferent authors. In what follows we summarize the main 
analytical results that can be found in the existing liter- 
ature. 

1) Tse et al.^ studied Coulomb drag between two 
unhybrydized graphene sheets separated by a dielectric 
by employing Boltzmann transport theory. They ne- 
glected the spatial dependence of the dielectric constant 
in the z direction (see Fig. [l]) and assumed a momentum- 
independent scattering time. In the weak-coupling limit, 
i.e. in the limit in which the interlayer distance d is 



much larger than the average distance between two elec- 
trons in each layer, Tse et al^ demonstrated that the 
low-temperature drag resistivity is given by 



Pd 



h ^C(3) (fceT)^ 



e2 32 eF,ieF,2 {qTF,id){qTF,2d) 
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{kF,id){kF^2d) 



(1) 



In Eq. (fl) gxF.f is the Thomas- Fermi screening wave vec- 
tor, which is proportional to k-p^e, the Fermi wave number 
in each layer. In the symmetric ni — n2 case the previous 
equation yields 



h T^ 



Pd oc 



^3^4 



(2) 



2) Wc now move on to summarize the main results 
by Peres et al.^'' . The authors of this work used Boltz- 
mann transport theory, took into account the momentum 
dependence of the scattering time and also the spatial 
dependence of the dielectric constant in the z direction. 
While the main approach followed by these authors is nu- 
merical, they do also provide an analytical expression for 
the Coulomb drag transresistivity in the weak-coupling 
regime. For an intralayer scattering time that depends 
linearly on momentum, Peres et alW^ found 



h T^ 



Pd fx 



l4d6' 



for fcpof ^ 1 



(3) 



3) KatsnelsorpsI used Boltzmann transport theory, 
took into account, at least partially, the spatial depen- 
dence of the dielectric constant in the z direction, but 
did not take into account the momentum dependence of 
the scattering time. The main results of Ref. HH! are: 



PD oc — -T^ J— ^ ^, for kpd < 1 

e^ n 

h T^ 
Pd oc -— -Ym, for kpd > 1 



(4) 



Note that the weak-coupling result by Katsnelson is in 
agreement with that of Tse et al^. 

4) Hwang et al^ used Boltzmann transport theory, 
took into account the momentum dependence of the scat- 
tering time but neglected the spatial dependence of the 
dielectric constant in the z direction. For a momentum- 
independent intralayer scattering time they find: 




h T^ 
h T^ 



, for k-pd <Si 1 
, for kpd ^ 1 



(5) 



Note that the weak-coupling result in Eq. (Is]) differs from 
the resuh by Tse et al0. 



For an intralayer scattering time that depends linearly 
on momentum, Hwang et al!^ find instead 



Pd oc -^ 






h T2 
Pd oc -^-j-ij, for kpd > 1 



(6) 



By comparing Eq. ([5| with Eq. ^ , Hwang et al^ con- 
cluded that the functional dependence of pD on n and 
d is very sensitive to the functional dependence of the 
intralayer scattering time on momentum. 

5) Narozhny et al^ have recently presented a sys- 
tematic study of Coulomb drag between two MD2DESs 
which is based on perturbation theory in the dimension- 
less coupling constant aoe = s^/i^v), v being the Dirac 
velocity. The authors of Ref. [29] did not consider the 
spatial dependence of the dielectric constant along the 
z direction and discussed mostly the case of "low dop- 
ing" [i.e. the regime in which the chemical potential is 
comparable with or smaller than k^T). Since the focus 
of our Article is on the opposite regime, i.e. the Fermi- 
liquid regime, we now provide a short summary of the 
results of Narozhny et alW^ in this regime only. The au- 
thors of Ref. [29l have demonstrated that the energy de- 
pendence of the relaxation time is completely irrelevant 
in relation with the low-temperature drag transresistiv- 
ity. In the weak-coupling limit Narozhny et al.^^ found 
for the drag transresistivity the same doping- and inter- 
layer separation-dependence as in Refs. [M] and [3S1 As 
density decreases, Narozhny et al^ found that the drag 
coefficient acquires logarithmic corrections - see Fig. 3 in 
Ref. [29l In particular, in the limit kpd <S^ 1, Eq. (41) in 
Ref. [2ni reads 



h T\ 
Pd oc - ^ — In 



Nfttcc 



(7) 



As we will see below, our results agree with those of 
Narozhny et al^, and represent their generalization 
to the case of a finite Thomas-Fermi screening length 
(oc A'facc) and spatially-dependent dielectric constants 
along the z direction. 



III. MODEL HAMILTONIAN AND BASIC 
DEFINITIONS 

We consider two unhybridized layers of massless Dirac 
fermions, each one described at the noninteracting level 
by the following single-channel Hamiltonian {h = 1): 



ni = vY^ V'^,^ (o-ap ■ k) i), 



k,l3,l 



(8) 



k,a,l3 



Here v is the bare electron velocity, k is the k ■ p mo- 
mentum, a, /3 are (sublattice) pseudospin labels, and 
Cap = i'^'aa^^aa) ^^ ^ vcctor of Pauli matrices which 
act on the sublattice pseudospin degree-of-freedom. The 



field operator Wka i i'4'k,a.i) creates (destroys) an elec- 
tron with momentum fc, pseudospin a, and layer index 
£= 1,2. 

In each layer, the Hamiltonian (Is]) can be easily diag- 
onalized by the matrix 



U{k) = 



1 



g-iVfc/2 g-JVfc/2 



(9) 



/9 \ e*'^'"/^ _g*¥'fc/2 

where Lpk is the polar angle of the vector k. 

The i-th Cartesian component of the current-density 
operator (z — x, y) in the £-th layer is given by 

JqA = Y. '^'i-q.aA''<p)^k,^^i ■ (10) 

k,a,l3 

The Fermi wave number in the £-th layer is defined by 



Ni, 



(11) 



where n^ > is the excess electron density in the f-th 
layep^S and Nf^i is a degeneracy factor {Nf^g = 4 for a 
graphene layer, accounting for spin and valley degenera- 
cies, while Nf£ = 1 for a TI surface state). 

In the absence of disorder, the Green's function corre- 
sponding to the Hamiltonian in Eq. ([8]) in the imaginary 
frequency axis is given by the following 2x2 matrix 

(/Zf -|- iuj)ta- + vk ■ cr 



Gi{k, iuj) 



(12) 



{m + iw)^ — w2fc2 

where 1^ is the 2x2 identity matrix in the sublattice- 
pseudospin representation and ^e is the chemical po- 
tential in the ^-th layer. In the zero-temperature limit 
fii — >■ Sfj = vk-p^i, where £f,£ is the Fermi energy in the 
f-th layer. 

The disorder-free Green's function in the eigenstate 
representation reads instead 

Ge,x{l^,ii^) == -. — , , (13) 

iu} + fie- ek,x 

where ek,x = Xv\k\ are Dirac-band energies. 

In what follows we will need the following matrix ele- 
ments: 

PkXM'X' = [uHk)U{k')]^^, 



(14) 



and 



<^kx,k'X' ^ [U\k)a^U{k')\^^, 

\Q-i{Vk + Vh')/'2 J^ ;\'g«(Vfc+¥'fc')/2 



(15) 



The two MD2DESS described by Eq. ^ are cou- 
pled electrostatically by long-range Coulomb interac- 
tions, which are influenced by the layered dielectric envi- 
ronment (see Fig.nl). The coupling Hamiltonian reads 



'^<=<= = ^Q XI ^U'{q)Pq. 



eP-q,e' : 



(16) 



where 



Pq,i 



EC 



q,cii,£ 



t 



k.a.- 



(17) 



is the density-operator for the £-th layer and Vu' (q) (with 
£ 7^ £') is the 2D Fourier transform of the interlayer 
Coulomb interaction 



^''^''^ = ^''^''^ ^ ^) '' 



(18) 



Here 



D{q) = [(ei+e2)(e2+e3)e''''+(ei-e2)(e2-e3)e""'] ■ (19) 

For future purposes, we introduce the dynamically 
screened interlayer interaction [/i2(g,i^), which, at the 
random phase approximation (RPA) level, is given 



where 



Ui2{q,(^) 



.(0) 



V^12(g) 

s{q,uj) 



(20) 



e{q,u) - [l-V,^{q)xriq,^)m~V22iq)xriq,^)] 



(0), 



~ V,\{q)x^°\q,i.)x^^\q,cu) 
is the RPA dynamical dielectric function. 



(21) 



In Eq. (21) 



^ {q, uj) is the well-knowrP^EH density- 
density (Lindhard) response function of a noninteracting 
MD2DEs at arbitrary doping ng. The Coulomb interac- 
tion in the £ = 1 (top) layer is given by 



Viiiq) 



qD{q) 



[(£2 + e^)e'i'' + (€2 - €3)6- 



-qd-\ 



(22) 



while the Coulomb interaction in the bottom layer, 
1^22 (?), can be simply obtained from Vii{q) by inter- 
changing £3 <->■ ei. 



Eqs. (18 1, (19), and (22) have first appeared in Ref. [ 



and their explicit derivation has been reported in Ref. 1261 
Notice that in the "uniform" ei = £2 = £3 = e limit 
we recover the familiar expressions Vii{q) = V22('?) — >■ 
27reV(eg) and ¥12(9) = V2i{q) -^ Vii{q)cyip{-qd). Most 
of the previous work on Coulomb drag in DLG has as- 
sumed this limit, which rarely applies experimentally. 

The aim of this Article is to present a theory of 
Coulomb drag, which is valid up to second order in the 
dynamically-screened interaction C/i2(g,w), for the sys- 
tem described by the Hamiltonian 



H-E^^ + ^° 



(23) 



q,t.^l' 



Note that we are not including in Eq. ( 23 ) any term de- 
scribing intralayer electron-electron interactions. These 
can be treated in an approximate fashion by invoking 
Landau's theory of normal Fermi liquidd^S]^ j_g_ ^y renor- 
malizing the microscopic parameters of the intralayer 



q, tuj„ + I 



fc+, ze„.+ + i^ 
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a) 



''fc,'i£„ k\'ie'.^H. 
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FIG. 2: Second-order Aslamazov-Larkin-type diagrams 
contributing to the Coulomb drag conductivity (finite- 
temperature Matsubara formalism). Solid lines denote the 
single-particle Green's function in the presence of disorder. 
Wavy lines denote the screened interlayer interaction U12 in 
Eq. ( |20[ |. Black dots in the triangular portions of the dia- 
grams denote vertices of current operators in the two layers, 
which are both proportional to the Pauli matrix a^ , say, if 
one is interested in the longitudinal drag conductivity. Fi- 
nally, k± — k ± q, fcj- = k' ± q, ie„.± — ie„ ± ia;„, and 



Hamiltonian i-Li. For example, in Eg. (|8J) one can use 
the renormalized quasiparticle velocity^^t'^, instead of 
the bare velocity v. Anyway, a treatment of the impact 



of intralayer interactions on the Coulomb drag transre- 
sistivity is well beyond the scope of this Article. 



IV. KUBO-FORMULA APPROACH TO THE 

CALCULATION OF THE DRAG 

CONDUCTIVITY 

Coulomb drag starts at second order in a perturba- 
tive expansion for the "drag conductivity" ctd in powers 
of the interlayer Coulomb interaction. To avoid infrared 
pathologies (stemming from the long-range nature of the 
Coulomb interaction) of certain integrals that appear in 
the theory. Coulomb interactions must be screened: from 
now on we will work with the dynamically-screened in- 
teraction C/12 introduced in Eq. ( 20 ) rather than with the 
bare potential V12 in Eq. (18 1. 



The two Aslamazov-Larkin-type diagrams^ that con- 
tribute to ctd up to second order in U12 are depicted in 
Fig. [2j In this figure, solid lines represent disordered 
propagators in the sublattice-pseudospin representation, 
while wavy lines represent the screened Coulomb inter- 
action Ui2. Finally, the vertices (black dots in the trian- 
gular diagrams) are current operators, one in each layer, 
evaluated at q = from the very beginning since we are 
interested in the drag conductivity in the uniform limit. 
Both current operators can be taken along the x direc- 
tion since we are interested in the current response of the 
"passive" layer in the x direction to an electric field ap- 
plied in the same direction in the "active" layer, i.e. we 
are interested in the longitudinal drag conductivity. 

Straightforward algebraic manipulations lead to the 
following compact expression for the sum of the two dia- 
grams in Fig. [2] 



0'D(«f^m) 



2^^ J W^^ 



y^ Ul2{q, iUJn)Ui2{q, iuJn + iftm)^l{q, i^n + if^m, iuJn)T2{q, iuJn, iuJn + i^r, 



(24) 



where j3 — {k-QT) ^ is the usual thermal factor, f7„j and w„ are bosonic Matsubara frequencies and 



Ti{q,iuji,iuj2) = ~^ X! / 77f^^^'[Giik,iSn)Gi{k + q,iSn + i(^2)<j"^Gi{k + q,ien + iuJi) 
+ Gi{k, ien)Gi{k - q, i£„ - ioji)cr''Gi{k - q, «£„ - 1^2)] 



(25) 



r 



is the so-called "non-linear susceptibility". Here e„ is a to the layer degree of freedom, and the symbol Tr[. . .] 
fermionic Matsubara frequency, the index ^=1,2 refers denotes a trace over sublattice-pseudospin indices. 




FIG. 3: The contour of integration C chosen to carry out 
the Matsubara sum in Eq. \27\ . Note the branch cuts on 
the real axis and in the lower half of the complex plane at 
Q'm(z) = —iQ.m- The crosses denote the poles of the Bose- 
Einstein thermal factor 71,3(2). 



We now need to p erfo rm the analytic continuation 
iVtm — >■ ri + zO+ in Eq. (24) and then take the d.c. fi — >■ 
hniit. It is very well known that the analytic continuation 



must be carried out after performing the sum over the 
Matsubara frequencies w„. Let us assume that Ui2{q, z), 
seen as a function of the complex frequency z, is analytic^ 
i.e. that it has no branch cuts. Under this assumption, 
the function T^{q^ z, z') has two branch cuts when z or z' 
are on the real axis and is analytic elsewhere, including 
the points z = and z' = 0. We then introduce 

/(q, z, z') = Ui2{q, z)Ui2{q, z'jTiiq, z' , z)T2{q, z, z') , 

(26) 
and perform the Matsubara sum in Ec^. ( 24 1 using a stan- 
dard procedure: 



Jit^r. 



- 15^fi9''^^n,i^r. 



- i^n) 



dz 



-.nBiz)f{q,z,in^ + z) , (27) 



2TTi 



where nB{z) = (e^^ — 1)^^ is the Bose-Einstein occupa- 
tion factor on the complex plane and C is an appropriate 
contour which excludes the branch cuts of the integrand 
(see Fig. Isl), located at 5to(z) = 0,—ifljn- Notice that 
the contour C as defined in Fig. [3] includes the points 
z — and z ~ —iQm, where the integrand is analytic. 
We obtain 



/+00 7 
---^nB{uj)\f{q,UJ^,ifl,n + Uj) - f{iVLrn,UJ^,i^rn + ^) 
-00 27rz 

+ /(if]m,w-zri™,a;+) - /(iri,„,a;- ir2„i,w")] , 



(28) 



since only the integrals around the branch cuts contribute 
to J{i^rn)^ while the integral over the circle vanishes 
when its radius is sent to infinity. Here uj^ — ijj ± iO^ 
and V stands for the Cauchy principal value, i.e. the 



point w = is excluded from the integration. 

We are now ready to perform the analytical continua- 
tion to real frequencies. After some algebra we find 



J{n + ie) ^vf -^\[nB{n + uj)-nB{uj)]f{q,uj-,n + uj+) 
J -00 27rz L L 

-I- nB(w)/(q',a;+,f} + a;+) - 713(1^ + w)/(q, w", rj + a;")| 



(29) 



We finally take the d.c. fimit fJ ^- in Eq. (29 1 and 
we obtain 



/+00 
-00 



dto 9nB(a;) 
27ri du) 



f{q,uj ,w+) , 
(30) 



since the terms in the last line of Eq. ( 29 ) vanish at least 
like ri^. Thus the drag conductivity in the d.c. limit 
read^ 



0-D 



IGtt 



(2^) 



+ CX3 



doj 



|£/l2(g,C^)P 

sinh2(/3w/2) 



ri(q,w^ 



)r2(q,w" 



(31) 



This is the central result of this Section and was first ob- 
tained in the context of Coulomb drag between ordinary 
2D parabolic-band electron gases - see e.g. Ref. El 



A. The non-linear susceptibility 

Let us go back to the definition of the non-linear sus- 
ceptibility r^ given in Eq. ( 25 ) . We introduce the follow- 
ing definition 



cedure: 



{q, le, le , le 



vN, 



dPk 



Tr[G,(fc,ie) 



(2^)2 
G,(fe + g,ieVG,(fc + q,ze")] 



(32) 



so that Eq. ( 25 ) reads 



= Jg{q,iuji,iu!2) 

+ Jg{-q,-ii^2,-i^i) ■ (33) 



We here remind the reader that uJi and W2 (^n) are 
bosonic (fermionic) Matsubara frequencies. Let us now 
concentrate on the first term on the r.h.s. of Eq. (33 1, 



i.e. on jTg{q, iwi, iuj2)- For the sake of simplicity, in what 
follows we omit to indicate explicitly the q-dcpcndcncc 
of the functions gi{q, is, le', ie") and Jg{q, iu!i,iuj2) every 
time we write equations that involve them. The complete 
notation will be restored only at the end of the mathe- 
matical manipulations that we perform on Eq. ( 33 1 . 



We first perform the fermionic Matsubara sum in the 



first term on the r.h.s. of Eq. (33 1 by the standard pro- 



Jg{ioJi,iuj2) = -'^gi{ien,i£n+iuJ2,i£n+iuJi) 



dz 
2Tri 



nF{z)ge{z, z + iu}2,z + iui) 



(34) 



where np^z) = {e^^ + 1)~^ is the Fermi-Dirac occupation 
factor on the complex plane. The function g(z, z+iui2, z+ 
iitJi) has three branch cuts in the complex plane located 
at 5m(z) = 0,—iu!i,—iuj2. Choosing a suitable contour 
of integration C, which encircles only the poles of nF{z) 
and leaves outside the branch cuts (in complete analogy 
with the contour drawn in Fig. ^, we find 



Jg{iuji,iuj2) 



+ 00 



de 



npie) [gi{e^,e + iw2, e + iuJi) 



, 27ri 

- gg{e',e + iuj2,e + ioJi) 

+ g^{e — iuJi,e^ + iuj2 ~ ii^i,£^) 

- gi{e — iuJi,£^ + iuj2~ ii^i^s^) 
+ ge{e - iuj2,e^,e^ +iu}i- iuj2) 

- ge{e - iuj2,e^,e^ +iuji- iuj2)] . (35) 

Here e^ = e±iO+ and we used that nF^s + icUm) = ^f(£) 
if LUjn is a bosonic Matsubara frequency. Notice that we 
are allowed to remove the index ± from e^ when it is 
summed to icoi or iu}2 (since they will be analytically 
continued to uj^) but not when it is summed to iuJi — iuj2. 



Analytically continuing icji -^ uj^ and iuj2 



[to 



obtain e.g. Ti{q,uj^ ,uj ) in the integrand on the r.h.s. 
of Eq. (31 1] we finally obtain 



^g('^l":^2 ] 



-t-oo 



de 
27ri 



np{e)[gi{e~^,e + uj2,s'^ + uji) - gi{e ,e +uj2,£^ +i^i) 



+ gi[e -wi,£ +uj2-uji,e'^) - ge,{e -wi,£ +W2-a;i,e ) 

+ giie'^ - W2,£^,£^ +^1 -t^2) - geis^ - W2,£",£" +UJi - W2)] 



(36) 



r 



We observe that the second term on the second line and with poles on the same half of the complex plane. After 
the first term on the third line of Eq. ( 36 ) are exactly zero some algebraic manipulations we get 



since they involve products of three Green's functions 



Jg{^t^^2 ) = 



+ C30 



de 

27ri 



{[nY{e) -nY{e + ijj2)]gi{e'^,e +a;2,e++wi) 



+ [n'p{e + uji) -nY{e)]gt{e ,£ + a;2,e+ + wi)} 



(37) 



We now take the limit uji,uj2 -^ lo and recast Eq. ( 34 1 in the following form: 



Jg{u}+,uj^) = vNf^e I —[nF{e + uj)-nF{e)] 



l-^Tr[[Gf{k,e)~Gf{k,e)]Gt{k + q,e + u;)a-Gfik + q,e + io)} 



(38) 



where we have introduced the retarded (advanced) treated in an analogous manner. The final expression for 
Green's function Gf^^\k, e). the non-linear susceptibility reads 



The second term on the right hand of Eq. ( 33 1 can be 






T,{q,uj+,LJ-) = vNi^eJ ^.[nF{e + io)-nF{e)]f-^Tr{[Gf{k,e)~Gf{k,e)]G^{k + q,e + Lj)a' 
X Gf(fc + q,£ + w)} + {(q,w)^(-q,-cj)} . 



(39) 



Eq. (39) together with Eq. (31) represent the most im- Coulomb drag transresistivity. 



portant results of this Section and are the starting point 
for the calculation of the drag transresistivity po ■ 



V. BOLTZMANN-TRANSPORT LIMIT 

Eqs. (31) and (39) need to be simplified for any prac- 
tical purpose and for a quantitative estimate of the 



We first switch from the off-diagonal sublattice- 
pseudospin representation of the Green's functions to 
the diagonal representation. In the latter representation 
Eq. (pgl) reads 



/+00 rip r rl^h f 

Y-- hr^ ^ {[Me + ^)-nF{e)][gtx{k.e)-gl^{k,e)] 

X gty{k + q,e + Uj)Qf,y, (fe + q, £ + w) Pk\,k+q\' ol+^y ^t,+q\" Pk+q\" ,k\j 

+ {(g,w)^ (-q,-w)| , 



(40) 



where we used the definitions in Eqs. (14)-(15) for the density and current vertices. 



We now show that the terms with A' ^ A" do not We then consider the term on the third hne of Eq. ( 40 1 



contribute to T£{q,uj~^,u! 

Let us in deed consider all the terms in the triple sum 
in Eq. (liol) in which A" = -A' = A'. We first recall that 



O" 



-, = iA'sin((/jfc+q) and 



k+qX' .k+qX 



PkX.k+qX' Pk+qX' ,kX 



iW 



sin((^fc - <Pk+q) 



(41) 



and perform the following changes of variables: k 
—k and e — > — e. Using the fact that sin((/3_fc_q) = 
-sin((/jfc+q) and sm{ip^k-ip_k-q) = sm{ipk - ipk+q) , we 
can write the off-diagonal (OD) contribution (A" = —A') 

torf^{q,u;+,u}-) as 



+ 00 



X [Gt\if^, e) - Gf,x{k, e)] G^x' {k + q,e + uj)Ql-y {k + q,e + uj) s\n{ipk - (pk+q) sm{(pk+q)} . (42) 



r 



Since n-p(z) + n-p{—z) = 1, the term in square brackets is 
identically zero. The OD contribution to r^(q,aj+,w~) 
thus vanishes. From now on we can set A" = A' in 



Eq. (40) and sum only over A and A'. 



We can further simplify Eq. (40) by assuming that, in 



the presence of weak disorder, the Green's function in the 
i-th layer can be well approximated by the expression, 



5Tik,u) 



^'"^^ 2T,(fc) 



(43) 



where Te(k) represents a momentum-dependent scatter- 
ing time in the i!-th layer and ^^ 3^ = Xv\k\ — pii are Dirac- 
band energies measured from the chemical potential fi£ 
of the ^-th layer. 

Note that here T£{k) should not be interpreted as a 
one-particle scattering time, but rather as the transport 
scattering tim^. This statement can be justified within 
the Kubo formalism by including disorder-related vertex 
corrections, or, in a much more transparent way, by using 
the Boltzmann transport equation, see Appendix [K\ 

Below, we will assume that the transport scattering 
time Ti{k) is isotropic, i.e. that it depends only on 
k — \k\, but allow different scattering times in the two 
layers. The specific functional dependence of ti on k de- 
pends on a particular model of intralayer impurity scat- 



tering. In the simple case of a momentum-independent 
scattering time we recover the usual "relaxation time ap- 
proximation" . A scattering time which depends linearly 
on momentum, T^(fc) oc fc, is a very popular model in 
the graphene literature. Early on it was understood^ 
that such functional dependence of Ti on k is needed to 
explain the linear-in-carrier-density d.c. conductivities 
that are experimentally measured in samples on dielec- 
tric substrates such as Si02 (and is typically attributed to 
charged impurities located close to the graphene sheet). 
Below, we will not assume any specific functional depen- 
dence of Ti on k. Our low-temperature analytical results 
for the Coulomb drag transresistivity do not depend on 
the particular scattering model that one chooses. 

At this point, it is useful to introduce the one-particle 



spectral function corresponding to Eq. ( 43 ) , 



Ae.x{k,u!) = ^mgf^{k,uj) 

n 

= t[gl^{k,uj)-gfjk,L,)] , (44) 

and the identity 

gf^{k,uj)g^^{k,uj) = n{k)Ai.x{k,oj) . (45) 



Using these definitions in Eq. ( 40 ) we obtain the following 



approximate expression for the non-linear susceptibility: 



Ti{q,uj+,uj ) w vN[^i 



;r~ / 7:^-T^}A[nFi£ + ^)-nF{e)]iTe{k + q)Ae.xik,e)At.x'{k + q,e + uj) 
-00 27ri7 {2TTy ^ L^ 



1 + AA' COs{ipk - ^k+q) X 



K + qX'.k + q,X'] + {{(l^^)^ (-9,-^)} 



(46) 



From now on we will introduce the simplified notation J^ ^ = va^-^ ^^.^ = vXcos(ipk), where in the last equality 
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we have used Eq. (15). The definition of J^ ^ should 
not be confused with tlie definition of the current- dens ity 



operator J^ ^ in second quantization given in Eq. (|10|) 



In this Article we are interested in the limit in which 
intralayer scattering is weak, which is, for example, the 
most relevant regime for high-quality DLG samples. In 
this limit we can approximate the spectral functions in 
Eq. (46) with 6 functions. 



At^x{k,£) w 5{uj- 






(47) 



Straightforward algebraic manipulations of Eq. (46) 



with the use of Eq. ( 47 ) yield the following Boltzmann- 



transport (BT) expression for the non-linear susceptibil- 
ity: 



Tf^iq,u;+,u;-) 



A, A' •' ^ ''^' 

X ...|-^g-^^^— ^[l + AAcos(,.,,-,.)]j 



(48) 



r 



We have dubbed Eq. (|48|) "BT expression" for the non- 



linear susceptibility since after inserting it in Eq. 
one obtains precisely an expression for the Coulomb drag 
conductivity that can be derived within a Boltzmann- 
equation approach. This will be shown in Appendix [A} 

Note that, in the usual "relaxation-time approxima- 
tion" (r^ independent of k), the BT expression for the 
non-linear susceptibility is simply proportional to the in- 
tralayer scattering time, r^"'"(<7,a;) ex r^. 



A. Coulomb drag transresistivity in the 
Boltzmann-transport limit 

We now focus our attention on the quantity which is 
actually measured in experiments, i.e. the drag tran- 
sresistivity pd, which is precisely the ratio between the 
voltage drop in the passive layer and the current in the 
active layer is indeed pD- This quantity can be easily 



found by inverting the 2x2 conductivity matrix, 

1 au 



PB 



-o-D- 



det 



CTi ctd 
CD 0'2 



crifT2 



(49) 



where a£ is the intralayer conductivity and the last ap- 
proximation in Eq. (49) holds true only if ctd ^ cf£. 



Within BT theory the intralayer conductivity at finite 
temperature is given bjB^l 



(7f ~ V Nf t 

2 



E 



(pk 

(2^ 



n{k) 






, (50) 



while the drag conductivity ctd is given by Eq. (|3lj) with 
the expression ( 48 1 for the non-linear susceptibility. The 



final expression for the BT drag resistivity reads 



BT 

Pd 



/3e2 



167rCTicr2 J (27r) 



-t-oo 



duj 



|t/i2(g,^)P 

■sinh^(/3a;/2) 



rr(q,^)rBT(q,^) 



(51) 



Note that in the relaxation-time approximation the BT 
drag transresistivity does not depend on the transport 
scattering times in the two layers (even if these are dif- 
ferent), since, as noted earlier in Sect.[V] F^""" is propor- 
tional to Tf in this approximation and so is ct^. 



The situation turns out to be much more complicated 
in the case in which r^ is an arbitrary function of fc = 
|A;|. A full numerical treatment of the three-dimensional 
integral in Eq. (51) at any finite temperature is beyond 



the scope of the present work and will be the subject of 
a forthcoming publication. 
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In the next Section, however, we demonstrate that a 
dramatic simphfication occurs in the low-temperature 
hmit. We will indeed prove that when k^T is much 
smaUer than the Fermi energies in the two layers, k^T <ti 
min£(eF,£), p%"^ is insensitive to the precise functional 
dependence ti — Ti[k) of the intralayer scattering time. 
Since for typical values of doping the Fermi energy in a 
graphene sheet is very large (w 1400 K for a carrier den- 
sity on the order of 10^^ cm~^), analytical results in the 
low-temperature limit are very useful when both layers 
are sufficiently away from the charge neutrality point. 

In the intriguing regime in which one layer (both lay- 
ers) lies (lie) at the charge neutrality point BT theory is 
inapplicable. Coulomb drag in this r egi me is dominated 
by unavoidable electron-hole puddles^ and is certainly 
very interesting but well beyond the scope of the present 
Article. For reasons of symmetry, lowest-order Boltz- 
mann transport theory applied to the regime in which 
one of the t wo layers is at the neutrality point gives 



VI. THE LOW-TEMPERATURE LIMIT 



In this Section we use Eq. ( 51 ) to calculate p^ analyt- 
ically in the low-temperature limit, for arbitrary values 



of the dielectric constants e.^ in Fig. [T] and for a generic 
scattering time Ti{k). 

As already anticipated above, the low-temperature 
regime is readily identified by the inequality k^T <C 
rah\i{eY^i)- In this limit we can: 

(i) use the low-t empe rature expression for the in- 
tralayer conductivitySSESl 



lim (Ji 

T-i-O 



47r 



Ni^ie'p,tTt{kF,t) ; 



(52) 



evaluate the non-linear BT susceptibilities 
(g, cj) and Ff "^(q, w) at zero temperature and only to 
lowest order in lj in the low-frequency u / min^ (ep.^) -^ 
limit; 



11 J 

pBT 
1^ 1 



and iii) replace the dynamically-screened interlayer in- 
teraction with the much simpler statically-screened in- 
teraction [/i2((7, 0). 

Indeed, the thermal factor sinh^^(/3cL'/2) in Eq. (51) 



represents an effective cut off on the values of uj in the 



integral in Eq. (51 1: it must be w ^ 2//?. 

To find the asymptotic behavior of F^'^(q,w) in the 



limit a; — > we first re- write Eq. ( 48 ) in the form 



rr(g,^) 



-TriVf 



/ Wf ^"'^^ ^ '^^'^'''''' " "'^''^•^''''^ f""^^^'^^^ " MCy)] 



X '^(^ + ^ki - d + q,A') [1 + ^A' COs(^fc + q - ifk)] 



(53) 



Clearly F^"'"(q,0) = 0. We thus need to calculate the 
derivative of rf'^{q,uj) with respect to cj at w = and 



in the zero-temperature limit. We find 



dTf^iq,u) 



doj 



1x1=0 



X.X'-' ^ ' 

X 5(e£ - Ca' ) [1 + AA' cos(^fe+q - ^fc)] 



T-i-O 



^ 5[np(^ + el,i)] 



duj 



'^" -nN,,e T.J^2Hk + g)4^,,, - n{k)Jl,] 5{i^^\)5{^':\ dt,A') 



A, A' 



X [1 + AA'cOs(^fc+q -(^fc)] , 



(54) 



where the last equality is valid only in the zero- temperature limit. More explicitly, we find 
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lim rF(q,a; 

T-i-O 



0) 



vr^iVw Y.I 4^ H^ + 9)'^^+'/.A - n{k)Jl^^]5{£,'i\)6{i 



fe+q,A 



') 



X [l + AA' cos{ifk+q - Vk)] 



(55) 



Eq. ( 55 ) is one of the most important results of this Sec- 
tion. ^An explicit expression for limT-j.orf"(q',a; — )• 0) 
will be given below in Eq. (58 1. Notice that the two 
^-functions in the above expression impose that only 
intraband excitations (A = A' = +1) contribute to 
T^^{q,uj — > 0) in the low-temperature limit. Moreover, 
the two (5-functions together pin the absolute values of 
k and fc -I- q to be equal to k-p.g- This implies that in 
the low-temperature limit the transport scattering times 
T(,{k + q) and Ti{k) inside the square bracket in Eq. ( 55 ) 
must be both evaluated on the Fermi surface of the ^-th 
layer, i.e. |fc -I- q| = |fc| = k^.i- Since Ti{k) depends only 
on the absolute value of its argument (and not on the 
polar angle of fc), both Te{k + q) and Ti{k) can he factor- 
ized out of the integral in Eq. \5^. T his is precisely the 
reason why at the end of Sect. |V A| we claimed that in 
the low-temperature limit p^ is completely insensitive 
to the precise intralayer scattering mechanism. 
More mathematically, we have: 

lim r^'^(<7,w — > 0) = --KLoNf iTilk-p e) 

T-)-0 



+00 



kdk 

~2^ 



6{vk - eF,i) 



dipk 



[Jk+q,+ Jk, + \ 



/o 27r 
X S{vk — v\k + q\) 
X [l + COs(<^fc+q - ipk)] . (56) 



We now employ the following identity, 

q 



JX 

•^k+q, 



_ JX _ 



C0s(v3q) , 



(57) 



which can be easily derived by using the aforementioned 
condition |fc -|- qj = |fe| = fcp/. The two integrals in 



Eq. ( 56 ) can be easily carried out analytically: we find 



lim Tf^iq, 

T-i-O 



0) = _^^ /'^y(fcF,^) e(2fcp, 

Zirvq 



q- 



K,i 



(58) 



where we have used that q^ — qcos{tpq). Eq. (58 1 is one 
of the most important results of this Section. We stress 
again that the applicability of this equation is not lim- 
ited to the case of a momentum-independent scattering 
time. Rather, Eq. (58 1 applies to a generic intralayer 
scattering time ti = Ti{k). The physical meaning of this 
result is that in the low-temperature limit the non-linear 
susceptibility is determined by what happens close to the 
Fermi surface. What matters thus is only the magnitude 
of Ti{k) evaluated at the Fermi momentum fcp/ of the l- 
th layer. This conclusion is in agreement with Narozhny 
et aim 



Using Eq. ( 58 1 and the following integral 



+ C30 



dx 



St:' 



smh'{yx/2) 3y^ 



(59) 



we finally arrive at the desired result for the low- 
temperature BT drag transresistivity: 



lim Pd 

T->0 



BT 



(fcBT)2 



6e^eF iEf 2V^ 



dq q\U,2{q,0)\' 



'1- 



4fc2 



'1- 



4fc|,2 



(60) 



where Qmax = niin(2fcF.i7 Sfcp 2). Once again, we stress 
that the result in Eq. ( 60 1 does not depend on the precise 
functional form of Ti{k). 



Dimensionless variables 



It is now useful to introduce dimensionless variables. 



We scale the wave number q in Eq. ( 60 1 with y^^kpskp^ 
by introducing x — g/-\/fcF,ifcF,2 and the effective in- 
teraction U12 with e"^/ ^ykpikp2 by introducing U12 = 
Ui2y/h\ik^/e''. 
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In these reduced units Eq. (60) reads (for physical rea- sons we restore Planck's constant from now on) 



r BT 

lini Pd = 



e^ 127r ep i£f 2 



dx x\Ui2{x, ( 



'1- 



x"^ k 



F,l 



4 fcp,: 



'1- 



c^ /c 



F,2 



4 /fc 



F,l 



(61) 



where 



mhi{2^JkY~iJk¥~2^ '^ \/ kp ,2 / kp ,1) ■ In is the most important result of this Section. 



Eq. (61) we have introduced the dimensionless coupling 
constant aoe = e^/(fi.w), which has a value ss 2.2 in DLG 
and « 4.4 in Bi2Te3 TIs if we use the respective Dirac 
velocities vq ~ 10^ m/s and v^i ~ 5 x 10^ m/s. Eq. (61) 



To proceed further, we write in an explicit manner 
the statically-screened dimensionless interlayer interac- 
tion i7i2(a;,0): 



Ul2{x,0) 



87ra;/i2(a;^) 



[X + 2iVi-iaeo\/fcF,l/^F,2 fllixm^ + 2A^f,2aeoVfcF,2/fcF,l f22ixO] " 16A^f ,liVt,2a?c/l'2 (^^O 



(62) 



where we have introduced the crucially-important dimen- 
sionless parameter 



^ = d^/kFJk^ , 
and the dimensionless form factors 

(e2 + £3)6^ + (e2 -£3)6" 



fiii^) 



/22(X) 



gix) 



(e2 + £i)e^ + (£2 "ei)e" 
.9(2;) 



/i2(2:) 



£2 

dix) 



(63) 



(64) 



(65) 



(66) 



with 



;(x) = (£i+£2)(£2 + £3)e'^ + (£l-£2)(£2-£3)e"'' • (67) 



The only ingredient we have used to write Eq. ( 62 ) is the 



well-known behavior of the static Lindhard function of a 
doped MD2DES for wave numbers q smaller than twice 
the Fermi momentum fcp ^: 



m 



xrii<2kF^i,o) = - 



Ni^ik-pj, 
2'irhv 



(68) 



Two asymptotic limits can now be easily inferred from 
the general low-temperature theory in Eqs. (61)-(67l 



1) the weak-coupling ("large interlayer distance" andTor 
"high density") limit, i.e. the limit in which ^ 3> 1, 
and 2) the strong-coupling ("small interlayer distance" 
and/or "low density" ) limit, i.e. the limit in which ^ <C 1. 



As discussed briefly above, the applicability of BT the- 
ory in the low-density regime is questionable: for the BT 
theory to remain valid, one should think of entering the 
strong-coupling regime by reducing the interlayer separa- 
tion d while keeping fixed the factor \/k-p,ik-p^2- Note also 
that, in the strong-coupling regime, diagrams of higher 
order with respect to those presented in Fig. [2] might be- 
come relevant. The applicability of RPA - Eq. ( pTj ) - is 
also questionable in this regime. Only a comparison with 
experimental results can tell us about the importance of 
these subtle issues of many-body theory. 



B. The weak-coupling limit 

In the limit i^ 3> 1, straightforward algebraic manipu- 
lations of Eqs. (61 )-(67) yield the following result for the 



low-temperature BT drag transresistivity: 



lim lim Pn 

{->-ooT->-0 



BT 



{k^Tf 



e" 3 ep^iep,2i^Nl^Nl^al, 

+00 



dy 



y'fUy) 



[hi{y)f22{y)-AfUyW 



(69) 



where we have introduced a new reduced variable, y — 
xS,. The quadrature with respect to the variable y can 
be carried out analytically yielding the following result: 



+00 



dy 



y'fhiy) 



[/ii(2/)/22(y)-4/2,(y)]^ 



+00 



dy 



sinh (y) 8 



c(3)£; 



(70) 
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where C{x) is the Riemann zeta function and C(3) w 1.2. 
Note that this quadrature does not depend neither on ei 
nor on 63. 

In summary, we find 



hm hmpgT = _A<(3) 
{->ooT->-0 e-^ 



8 d^alNl,Nl,kl,kl, 



£F,ieF,2 



oc — ^ £9 



(fcBr)2 



3/2^3/2 



d4 



(71) 

The dependence of /Oq""" on temperature, densities and 
intcrlayer distance displayed by Eq. ( 71 1 is in agree- 
ment with the results by Tse et al^^ and by Katsnel- 
son^SI, even though the numerical prefactor we find is 
four times smaller (we remind the reader that for DLG 
^f.i — Ni,2 = 4). It is important to observe that the 
weak-coupling and low-temperature BT Coulomb drag 
transresistivity is sensitive only to the dielectric constant 
£2 between the two MD2DESs (see Fig.[l]). 



C. The strong-coupling limit 



In the limit ^ — >■ we can expand all the functions 



fij{x£,) that appear in the effective interaction (62) is 



powers of their argument since x is bounded from above, 
< a; < Xmax, and ^ — ?> 0. Following this procedure we 
find the following asymptotic expression for [/i2(a;,0) in 
the hmit ^ — > 0: 

An 
limt/i2(a;,0) = 



{ei+e3,)x + 2 



^TFa + 9TF,2 



(72) 



yjk-psk 



F,2 



where we have introduced the Thomas-Fermi screen- 
ing wave number qj^i = -^f/Qioe^F/- Using this re- 
sult in Eq. (61 1 we find the final expression for the 



low-temperature BT drag transresistivity in the strong- 
coupling limit: 



lim lim pp, = 

5^0 T^O 



h An 



' £f.i£f.2 




X^ fcF,2 

4 fcp 1 



(ei-|-e3)2; + 2 



^TF,! + '7TF,2 



y/kF,lk 



F,2 



(73) 



Eq. ( 73 ) simplifies considerably if we assume equal den- 
sities and degeneracies in the two layers: ni — n2 = n, 



i-Fa 



«F,2 



^F, £f,i = eF,2 = £f, and A^^f^i 



N, 



f,2 



iVf. In this case Eq. (73|) reduces to 
h A-K 



lim lim Op, 



BT 



fcer 



F{ei+ez.,Nfac 



(74) 



where 



F{ei -|-e3,7Vfaoo) = / 
Jo 



dx 



a;(l-a;V4) 



[(ei + ez)x + 47Vfaee] 

(75) 
Note that in the strong-couphng hmit hmr^o pW' does 
not depend on the interlayer distance d and scales like 
1/n. As far as the dielectric constants e^ are concerned, 
the strong-coupling and low-temperature BT Coulomb 
drag transresistivity depends only on ei -f 63 but not on 

£2- 



The quadrature in Eq. ( 75 1 can be easily carried out 
analytically. The result is 



F{ei +e3,iVfaee) 



1 



2(ei + e3)^ 



12A^facc(ei+e3) 



- (ei + 63)2 [3 + 2 ln(2)] + 2ANfa, 
X ln(2) + 2[(6i+e3)'-12iV2a2j 

ei + £3 \ 1 



X In 2 



Nfa^ 



(76) 



Note that T{ei -I- e3,A'facc) diverges logarithmically in 
the weak-screening Njacc -^ limit, in agreement with 
Eq. (41) of Ref. M 



VII. DEVIATIONS FROM THE 

LOW-TEMPERATURE QUADRATIC BEHAVIOR 

DUE TO A FREQUENCY-DEPENDENT 

SUBSTRATE DIELECTRIC CONSTANT 

In this Section we illustrate how deviations from the 
Fermi-liquid quadratic-in-temperature dependence of the 
Coulomb drag transresistivity can occur in the situation 
in which the dielectric constant e^ of one of the media 
surrounding the two Dirac-fermion layers has a strong 
frequency dependence. 

Consider, for example, a DLG system on a substrate 
like SrTi03 . Following Ref. |45l we model the dielectric 
constant of this substrate by a frequency-dependent func- 
tion of the form 



f^si^) 



+ (eo - Coo) 



■ 170; 



(77) 



At room temperature Coo = 5.2 , eg — 310, ujo/{2tt) = 
2.7 THz, and 7/(27r) = 1.3 THz. 

We would like to understand what is the qualitative 
role played by the change €3 — >■ e3{uj) in the Coulomb 
drag transresistivity. For simplicity, we assume that the 
two graphene sheets comprising the double layer have 
the same density ni — n2 — n (/cp,! — kp2 = ^f)- We 
also assume that the Fermi energy ep corresponding to 
n is the largest energy scale in the problem: under this 
assumption, we are still allowed to expand the BT non- 
linear susceptibility Tf"^{q, lu) to lowest order in the small 
parameter huj/ep. 
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Despite Tiuj ^ ep, we can have two distinct regimes 
since £3(0;) brings in a new frequency scale, i.e. ujq: 1) 
w ^ Wq ^ £f/^ and ii) wq ^ t^; <C ep/h. In the first 
regime what matters is obviously £3(0; = 0) — eg. In 
the second regime, instead, we have to retain the full fre- 
quency dependence of £3 (w) . The existence of this second 
regime ensures the possibility to observe deviations from 
Pp""" ex T^ above a certain temperature scale. 

The low-temperature Coulomb drag transresistivity for 
the situation described above is given by Eq^51) with 
rf^ {q,uj) given by the expression in Eq. (58 1 and ag 
given by the expression in Eq. ( 52 ) : 



lim pj3 



BT 



h/3 



e^ Wn^v'^ep 



2fcF 



dq q[l 



4A;| 



+- , iJ^\U,2iq,0)\ 



doj 



£3-i-£3(w) 



sinh^(;3w/2) 



.(78) 



Here the notation Ui2{q, 0)|^ ^^ ,^-. means that from the 
point of view of the electronic system we still have to use 
the statically-screened interlayer interaction, while from 
the point of view of the substrate we have to take into 
account the frequency dependence of 63 through the use 



of Eq. (771 



The double integral in Eq. (78) can be easily per- 



formed numerically and some illustrative results are sum- 
marized in Fig. PI For temperatures T ^ huj^/k^ we 
find the usual Fermi-liquid quadratic-in-temperature be- 
havior, —p^/po = aT^ where po = h/e^ and a is a 
numerical coefficient whose actual value depends on the 
carrier density n and on the interlayer distance d. How- 
ever, we clearly see from Fig. [4] that, for temperatures 
T ^ huJo/kB, the low-temperature Coulomb drag tran- 
sresistivity deviates from the canonical Fermi-liquid be- 
havior. 

Further deviations from the Fermi-liquid temperature 
dependence are induced by the explicit (and strong) tem- 
perature dependence of the dielectric constant of SrTi03 
- see, for example, Ref. [31] - which here has been ne- 
glected for simplicity. 



VIII. SUMMARY OF OUR MAIN RESULTS 
AND CONCLUSIONS 

In summary, we have presented in a complete and as 
much pedagogical as possible manner the minimal theory 
of Coulomb drag between two spatially-separated two- 
dimensional (2D) systems of massless Dirac fermions, 
which are both away from the charge- neutrality point. 
Our theory relies on second-order perturbation theory 
in the screened interlayer interaction and on Boltzmann 
transport theory. 

In this well-established theoretical framework, which 
has also been adopted by other authors in earlier 
workPfflSSl, we have clearly demonstrated that the pre- 
cise functional dependence of the intralayer scattering 
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FIG. 4: (Color online) Top panel: the Coulomb drag tran- 
sresistivity (y^ (in units of — po = —h/e^) is plotted as a 
function of temperature T (in units of the Fermi temper- 
ature Tp = ep/fca). Note that both axes are in logarith- 
mic scale. Data in this panel refer to: a carrier density 
n = 3.1 X 10^^ cm^^, an interlayer distance d ~ 2 nm 
(dfcp = 5.0), £1 = 1, and £2 = 7.8. Filled circles label the data 
obtained by taking into account the frequency dependence 
£3 = £3(1^) reported in Eq. (77 1. The solid line labels the 
results obtained by setting £3 = £3(0; = 0) = £0 in Eq. (|78|. 



Bottom panel. The Coulomb drag transresistivity pp (in 
units of — poffl, where the coefficient a changes with changing 
dfcp - see text) is plotted as a function of T/Ty. Note that 
also in this panel both axes are in logarithmic scale. Curves 
labeled by different symbols correspond to different values of 
the interlayer distance d (density is fixed at the same value 



used in the top panel, i.e. n = 3.1 x 10^'' cm ^ 



time on momentum plays absolutely no role in determin- 
ing the low-temperature Coulomb drag transresistivity 
in the Fermi-liquid regime. This is in contradiction with 
the findings reported in Refs. [17] and [55] while it agrees 
with the conclusions reached by Narozhny et alW^. 
For two layers with identical degeneracies (iVf^i = 
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Ni^2 = -^f) and densities (ni = n2 = n), we find the 
following results for the low-temperature Coulomb drag 
transresistivity: 

1) In the weak-coupling limit - Eq. (71 1 -we find 



lim Pd 



BT 



h 7rC(3) 



k^T^ 



Sal^N^kl n\^ 



(79) 



The dependence of this result on layer-separation d and 
doping n agrees with that found in Refs. 12412^ and E51 
2) In the strong-coupling limit - Eq. (74) - we find 



lim Pd 

T-i-O 



BT 



e2 3 fc| 






-e3,iVfaeo) , (80) 



where the explicit expression of the function J^(ei -|- 
e3,N{aee) is reported in Eq. (76). The independence of 
this result on layer-separation d and the dependence on 
doping n (oc 1/n) agree with the findings of Ref. |29l 
even though the authors of this work captured only the 
weak-screening A'facc — >■ asymptotic behavior of the 
function J-{ei + £3, A'faoo) and did not take into account 
the spatial dependence of the dielectric constant along 
the z direction. 

General results for ni 7^ 712 and iVfi ^ N[,2 can be 
found in Eqs. (FtTI) and (l74|). 



Finally, we have shown how deviations from the Fermi- 
liquid temperature dependence can occur when one of 
dielectric constants of the media surrounding the Dirac- 
fermion layers depends strongly on frequency and tem- 
perature. Double-layer graphene systems deposited, for 
example, on SrTiOa , a v ery we ll-known insulator close to 
a ferroelectric instabilitjE^ES, can be used as a testbed 
for this idea. 

In the future, we plan to present a systematic numeri- 



cal study of Eq. (51 ) at arbitrary temperatures and to in- 



vestigate more deeply the strong-coupling limit by study- 
ing i) the impact of diagrams of order higher than two in 
the screened interlayer interaction and ii) beyond-RPA 
corrections. 
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Appendix A: Boltzmann transport approach to the 
drag conductivity 

The Boltzmann transport equation for the distribution 
function fe{k, A) of the £-th layer reads: 

dtfiik,\)+Vkfiik,X)-dtk = I,{k,X)+jfjik,X) , (Al) 

where Ie{k,X) is the so-called collision integral, which 
describes intralayer scattering events, while the term 
7^ g{k, A) describes momentum transfer between the two 
Dirac-fermion layers {(. = 2 \i i — 1 and £ — 1 ii £ = 2) 
and is obviously absolutely crucial in the Coulomb drag 
problem. We are interested in the steady-state regime 
in which dtfi{k,X) = and employ a "generalized re- 
laxation time approximation" for the intralayer collision 
integral, i.e. 



Ii{k,X) = - 



fi{k,\)-ff\k,X) 



Te{k) 



(A2) 



where f^ ' (fc. A) 



nA& 



is the equilibrium distribu- 



tion function. Below we will assume that the transport 
scattering time Ti{k) is isotropic. 

The interlayer collision integral reads 



d^k' 



,,,ik,X) ^ 2.N,, Y. /7|^/Sl^-(^'^i'l-C.v)P<l-CS.v + 
A',A",A"'-^ y > J y ) 

X {/,(fc, A)[l - Mk + q, X'Mik', A")[l - Mk' - q, A'")] 
- [l-/,(fc,A)]/,(fe + q,A')[l-/,-(fc',A")]/,-(fc'-g,A"')} 






^k'-q,X"') 



1 + AA' COs((^fe - (fik+q) 1 + A" A'" COs(l^fc' - Lpk'-q) 



(A3) 



r 



The degeneracy factor Nf i can be understood from a In the semiclassical limit dtk — —eEg, where Ei is the 

careful analysis of the conservation of spin and valley electric field in the ^-th layer. In the regime of small 
degrees of freedom during the scattering process. electric fields we can linearize the Boltzmann transport 



17 



equation by writing 



(A4) 



The linearized Boltzmann equation in the steady-state 
regime reads 



e\7kf^°\k,X)-Ee 



/W(fc,A) ^ ninca.)^^^^^^ 



Tiik) 



(A5) 



where the hnearized interlayer colhsion integral reads 



7, 



(linear) 



\i \ri Mil '^ \ / ^ 






til) 



tie) 



A', A". A'" 



[u 



+ f, 
- /, 



p(l 



p(l 



+ / 



(1 



(27r)2 7 (27r)2'"''^^'^'^fe.-^ '^*=+'J 

^) (fc, A)/f' (fc', A") [1 - /f ' (fc + q, A')] [1 - /f (fc' - q, A 
(fc + q, A')/f^ (fc, A)/f ) (fc', A") [1 - /f ' (fc' - g, A'")] 



?fc+<J,A'il "^?fe,A ';fe+q,A' ^"Sfc',A" ?fe'-q,A"V 



i"A-f(0)^ 



(0)/ 



(fc', A")/,^^'^ (fc, A) [1 - /,^"^ (fe + q, A')] [1 - /f ^ (fc' - q, A'")] 
(fc' - q, A"')/r(fe, A)/f (fe', A")[l - /r(fe + q, A')] 
(fc, A)/f (fc + q, X')ff\k' - g, A"')[l - /f (fe', A")] 
(fe + q,A')/f(fc'-q,A"')[l-/r(fc,A)][l-/f(fc',A")] 

(fc', A")/l°^ (fc + g, A')/f ^ (fc' - g, A'") [l - /f^ (fc, A)] 

(fe' - g, A"')/r(fc + g. A')[l - /^(fc, A)][l - /f (fc', A")]} 



1 + AA' COs{ipk - •fk+q) 1 + A" A'" COs(l^fc' - ipk'-q) 



(A6) 



We now seek solutions of Eqs. (A5)-(A6| of the form 



and 



p(i) 



(intra) / 



+ ^f "■)(fc,A)fc•£;,- 



(A7) 



_^^f-)(fc,A) , S,^f-\k,X) 



Ti{k) 



5Ef 



(A9) 



In the following we will assume that \'ipi^^^ {k,\)\ ^ 

\'ij)f^^'''^ \k^\)\. The coefficients V')'" "^ (fc, A) and 

V'f " ""^ (^i A) can be found by calculating the functional 
derivative of the linearized Boltzmann transport equation 
with respect to E^ and Eg, respectively, i.e. 



In writing Eq. ( A8 ) we have neglected the small term 
hfT'"'\kA)/5Ei. The solution of Eq. JasI reads 



eVfe/f^(fe,A)~-fe 



4'"^"^(fc,A) 
n{k) 



^r->{k,X) = ek ■ Vki^k\ ^^'t'^^ n{k) . (AlO) 



, (intra) , 



The second term on the r.h.s. of Eq. (A9| can be cal 



(A8) culated from Eq. (A6). Solving Eq. (A9) we find 



V', 



(inter) 



{k,\) 



\i \ii Mil J \ J J 



<fik' 



-\u..{,,i\l\-i\l' -' 



A', A", A" 



'^(eS - Cq^' + &' - eil,.^^0{/r (fc, A)/f (fc', A")[l - /r (fc + 9, A')][l - /f (fc' - g, A'")]} 



1 + AA' cos{ipk - Vk+q) 1 + A"A"' cos{ipk' - tpk'-q) 



To obtain Eq. (All) we used the identities 
9/f)(fc,A) 



e(o)i 



f(") 






/3/r(fc,A)[l-/r(fe,A)] (A12) 



and 



(A13) 

The current density in the ^-th layer reads 



Ji - -eiV,,5:|A^/«(fc,A)V.eS 



J (intra) j'(intcr) 



(AM) 



Using Eq. (|A10|) we find J^ 

. ^ (2.) 

f(0) 



(intra) _ 






• Ep 



d^i'i 



v.d^l 



e^A^f^f 



E 



d^fe 



2 ^7 (2^)2 



f(0) 



T£(fc) 



'<^ 
afc 



g/r(fc.A) 



(A15) 



In the last equality of Eq. (|A15|) we used the fact that 

(A16) 






dk 
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(All) 



(since S,k-\ depends only on |A;|) and that the average 
over the angle of (fc • v)k = v/2 for any constant vector v 
(in D = 2 spatial dimensions). From Eq. (A15) we find 



immediately the longitudinal intralayer conductivity 



5J, 



{i) 



a-£ 



5E^ 



(*) 






(Pk 



rdk)\ 



dk ] 









d!f\k,X) 



dCk 



w 



,(A17) 
where ? is a Cartesian index and in the last equali ty we 
used the fact that £,l^\ = Aw|fc| - /i^ Note that Eq. ( A17) 



coincides with Eq. ( 50 ) in the main text 



The interlayer current density reads 



jf-) = _eiv,,^| A^v('""^'(fc, A)(fc.£;,)v.ei'i , 

. . (A18) 

where ■0'^'"*'"'' (fc, A) is given in Eq. (All). After some 
straightforward algebra we obtain 



T(inter) _ „ 2nNf^iNf^^(3e 
'i e — I^e 7. 



E 



A, A', A", A" 



d^q f d^k f d^k' ,^^ ,_ ^(f) (£) ^|2 



i2ny J i2ny J (2^)2 '"^^^l?, 5fc,A Wg.A'^i 

SiCi'i - ^i'L. + d^lA" - dl,a'"){/r (fc, A)/f (fc', A")[l - /r (fe + 9, A')][l - /f (fc' - q, A'")]} 
ri{k)Wk^i'\-n{k + q)Wk^l^: 



q,\' 



,{£) 



ri{k')Vk'ek'[y' - ri{k' - q)Vfe,Cr'-^,A"' 



W 



1 + AA' COs((pfc - (/3fc+q) 1 + A" A'" COs{ipk' - fk'-q) 



(A19) 



Notice that Eq. (A19) has been recasted in a symmetric form by means of a simultaneous exchange of dummy 
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variables {k -^ k + q, k' ^ k' — q, \ -^ \' , and A" O A'") 
and thus the integral is multiplied by an additional factor 



1/2. From Eq. (A19) we can derive the longitudinal drag 



conductivity ctd, which reads as following: 



5J, 



(*) 



CD 



5E 



(i) 



16n 



(2^)2 J_ 



+ CX3 



dcj 



|t/i2(g,^)P 

■sinh^(/3a;/2) 



r-^(q,^)r-^(q,c.), 



(A20) 



where the function r^"''(<7,a;) coincides with the o ne de- 
fined in Eq. (48) of the main text. In writing Eq. ( A20) 



we have introduced an auxiliary variable w to disentangle 
fc from k': 



S{i - e') = / dcjSi^ - u:)S{e - io) 



(A21) 
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